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The organization in brain networks shows highly modular features with weak inter-modular inter¬ 
action. The topology of the networks involves emergence of modules and sub-modules at different 
levels of constitution governed by fractal laws. The modular organization, in terms of modular 
mass, inter-modular, and intra-modular interaction, also obeys fractal nature. The parameters 
which characterize topological properties of brain networks follow one parameter scaling theory in 
all levels of network structure which reveals the self-similar rules governing the network structure. 

The calculated fractal dimensions of brain networks of different species are found to decrease when 
one goes from lower to higher level species which implicates the more ordered and self-organized 
topography at higher level species. The sparsely distributed hubs in brain networks may be most 
influencing nodes but their absence may not cause network breakdown, and centrality parameters 
characterizing them also follow one parameter scaling law indicating self-similar roles of these hubs 
at different levels of organization in brain networks. 


One of the most important issues in the study of brain 
networks is the origin of functional modules, organiza¬ 
tion of these modules and their functional relationships. 
Brain networks, constructed from various experimental 
studies on brains of different species, exhibit hierarchical 
features (highly modular structure) [1, 2], and these mod¬ 
ules are believed to be sufficiently isolated to enable them 
to perform independent functions [3] . These sparsely dis¬ 
tributed modules in brain network are shown to exhibit 
small-world topology, which have large local clustering 
co-efficients and very small path lengths [4], and it may 
allow the modules to perform independent functions [5]. 
On the other hand, studies on brain networks derived 
from functional magnetic resonance imaging (fMRI) [6- 
8], structural MRI [9] and diffusion tensor MRI [10] show 
small worldness in brain networks which seems inconsis¬ 
tent with the observed high modularity. High cluster¬ 
ing in small worldness, which is a local parameter, could 
not explain the global high modularity of brain networks 
[1, 11, 14]; and short path length in small worldness is 
also not suitable for strong modularity [12, 14]. Since 
the strong modularity corresponds to large world, the hi¬ 
erarchically organized, highly clustered, nearly isolated 
and self-similar set of modules are shown to be weakly 
tied among themselves [12], as a consequence of which 
the network preserves the small-world properties [14]. 
Therefore, the weak ties among the modules are believed 
to maximize the information transfer among the mod¬ 
ules with minimum wiring cost, and also allow to main¬ 
tain small-world topological characteristics [14, 15]. Fur¬ 
ther, these weak connections among modules in brain 


‘Electronic address: brojen@jnu.ac.in (Corresponding author) 


networks compel limited propagation of avalanche of neu¬ 
rons among the modules and are modular size dependent 
[16]. 

Fractality or self-similar structures in a complex net¬ 
work could be one property which can explain functional 
relationships of a larger network down to the funda¬ 
mental structure through different levels of organization 
[12, 13] . Scaling and renormalization theory can probably 
highlight the importance of information flow in a complex 
network and its self-organization [17]. It has been shown 
that hierarchical organization of modules in functional 
brain networks (fMRI) show fractal properties [14]. This 
fractal organization of modules in the network keeps hubs 
tightly bound inside corresponding modules, and use low- 
degree nodes as inter-modular connectors showing disas- 
sortative topology [18]. However, whether organization 
of modules and sub-modules at different topological lev¬ 
els follow fractal nature or not is still an open question. 
Further, whether the fractal properties exist in brain net¬ 
works of different organisms (lower to higher level organ¬ 
isms) is not fully investigated. 

In this study, brain networks of three different species 
(from lower to higher level of brain organization), namely, 
C. elegans, cat, and macaque monkey were investi¬ 
gated. Scaling laws and fractal rules were applied on 
several topological parameters to investigate the self¬ 
organization and fractal properties of the brain networks. 

Results 

The topological properties of hierarchical network, which 
involve emergence of well-defined modules with few 
sparsely distributed hubs, can be characterized by three 
topological parameters, namely, probability of degree dis¬ 
tribution P{k), P{k) ^ k~'*, with 7 < 2.0 [19], clus- 
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FIG. 1: Hierarchical organization of brain networks of C. elegans, cat, and monkey at different levels. The upper parts show 
the topological arrangement of modules and sub-modules at various levels of organization (one way of largest module and 
sub-module) till the motif level. The lower parts show the organization of all modules and sub-modules at different levels 
(levels are indicated by circles) of the brain networks of the three species. 


tering co-efficient C'(fc), C{k) ~ with a Ri 1 [20], 
and neighborhood connectivity Cn{k), Cn{k) ^ k~^ with 
/3 Ri 0.5 [21]; and follow power-law distributions with de¬ 
gree k [22]. The set of the exponents ( 7 , a, /3) of the three 
distributions (T’(fc), C{k), C'„(fc)) calculated using net¬ 
work theory for the three species, namely, C. elegans, cat, 
and monkey are given by: C. elegans (2.0,0.65,0.28); 
cat ^ (1.8,0.8,0.25); and monkey —>■ (1.7,0.73,0.2), 
showing hierarchical features in all the three brain net¬ 


works studied (Fig. 1 and Fig. 2 panels (A)). The or¬ 
ganization of modules and smaller modules at different 
levels {level-1 modules are the set of modules constructed 
from network, level-2 sub-modules are the set of all sub- 
modules constructed from level-1 modules, and so on) of 
the three species (Fig. 1) shows the hierarchical prop¬ 
erties. The smallest module in each brain network (C. 
elegans, cat, and monkey), from which all the three dis¬ 
tributions {P{k), C{k), and Cn{k)) can be calculated, is 
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found to be triangle motif (Fig. 1). So one can think 
of triangular motif as the fundamental regulator of each 
brain network. 

Scaling of modules at different levels 

We now study the topological properties of modules at 
different levels in each brain network (Fig. 2 panels (A) 
and (B)). The plotted P{k) versus k of C. elegans of 
larger modules and sub-modules at different levels (Fig. 
1: level-1, C3; level-2, SC13; level-3, SC13-3) show nearly 
parallel straight lines of power-law fits of P{k) with k at 
log-log scale. The data points of the modules and sub- 
modules at different levels along with data of whole net¬ 
work are rescaled to a single plot using one parameter 
scaling theory (see Methods) [23-25] and the power-law 
fitting of P{k) on the rescaled data gives 7 = 1.8. This 
scaling behavior is satisfied to all the modules and sub- 
modules in the C. elegans brain network, but shown only 
for one path as shown by arrows in Fig. 1. The same pro¬ 
cess is done to calculate clustering co-efficient C{k) and 
neighborhood connectivity Cn{k) of the same data set of 
modules and sub-modules at different levels, rescaled the 
data set with original whole network using this one pa¬ 
rameter scaling method, fitted with power-law equations 
with k in log-log scale and found their exponents to be 
a = 0.27 and /3 = 0.7, respectively. 

Similar network is done for cat and monkey brain net¬ 
works also (Fig. 2). Scaling of the set of data of each 
species is done using the one parameter scaling method, 
and exponents of the respective distributions which spec¬ 
ify topological characteristics of the networks are deter¬ 
mined by fitting the rescaled data with the respective 
distribution equations (Fig. 2). The set of exponents 
of the distributions of the scaled data of the two species 
are: cat —>■ (1.7, 0.7, 0.25) and monkey —>■ (1.6,0.68,0.25), 
respectively. 

Fractal nature of modules at different levels 

The characterization of self-similar structures in network 
can be studied from the evolution of structures (num¬ 
ber of nodes in the structures) in the network with path 
length [12]. We calculated the number of nodes n and 
diameter i?L in each module or sub-module in a certain 
level L, and then average over number of nodes N = (n) 
and path lengths = (Rl) of all modules are taken. 
The evolution of N{r]^) as a function of for all levels 
in each network of the three species is shown in Fig. 3. 
The behavior of iV(rL) with in all the brain networks 
follows the following power law: 

( 1 ) 

where i = 1,2, 3 and d = {di-,i = 1,2,3} are Hausdorff 
fractal dimensions of brain networks of the three species. 
The value of di of brain network of each species can be 
calculated by fitting the power-law equation ( 1 ) on NI^vl) 
versus data of the respective species (Fig. 3), and the 
calculated fractal dimensions of all the three species are 
given in Fig. 3. The fractal dimension is found to be 


the largest for C. elegans {dee = 3.47) and smallest for 
monkey {d^ = 1.93). Since fractal dimension is directly 
related to surface morphology of any system, larger value 
of fractal dimension may probably indicate larger disor¬ 
der in network organization [26]. Its smaller value in 
the brain network of higher level species may reveal the 
organization of the network is more ordered and system¬ 
atically self-organized [27]. 

To understand our claim of fractal nature of organi¬ 
zation of modules (relating to the interaction) in brain 
networks, we now calculate the number of edges (e) and 
diameter {tl) in each module of a certain level (for jth 
module of Lth level: , r^^), and then obtain average 

edges and diameter of the modules of the level L of ith 
species given by A, = A and vl = ^ EJ" ’’l > 

where m is the number of modules/sub-modules at level 
L. The behavior of Ei as a function of again obeys 
the following power law (Fig. 4): 

E,{rL)~rf£, ( 2 ) 

where fi = {fi',i = 1,2,3} is the set of fractal dimension 
relating to edges of modules and sub-modules of brain 
networks of the three species. The fractal dimension val¬ 
ues of the respective species in this case are found to be 
higher than the respective values fractal dimension cal¬ 
culated using network mass or network node number, i.e. 
fi > di] however, both di and fi show the similar nature 
(Fig. 3 and Fig. 4 lower right panels). 

The nature of organization of modules among dif¬ 
ferent levels can also be investigated by studying the 
inter-modular interaction among the modules and sub- 
modules. We calculate the number of edges between any 
pair of modules in a particular level L of brain network of 
ith species, average over all the inter-modular edges of all 
possible pairs of modules/sub-modules given by F^, and 
then study the variation of Fi as a function of average di¬ 
ameter of all modules/sub-modules in the level (Fig. 
5). The variation of F^ with for all brain networks of 
the three species C. elegans, cat, and monkey (Fig. 5) 
shows power-law behavior (fitted line to the data points) 
given by 

r^{rL)-r|:, (3) 

where gi = {gi',i = 1,2,3} is the set of fractal dimen¬ 
sions of brain networks of all the three species. This 
power-law nature reveals the fractal nature of the inter¬ 
modular organization of the brain networks. The power- 
law behavior of mass (number of nodes), intra-modular 
and inter-modular edges of modules and sub-modules in 
all the levels of brain networks show the fractal organi¬ 
zation of brain networks. 

The scaling and fractal properties of modules and sub- 
modules at different levels of the brain network of each 
species probably connect the topological organization of 
the modules and sub-modules to their functionalities and 
working relationships among them, within and among the 
levels. Further, self-organization among these modules 
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FIG. 2: Topological characteristics of brain networks of three species, C. elegans, cat, and monkey: (A) for whole brain network, 
(B) for modnles and sub-modules at various levels of network, (C) scaled for all modules and sub-modules in all levels to a 
single plot. The first three upper rows of panels are for clustering co-efficient, next three rows are for connectivity, and last 
three rows are for probability of degree distribution of the three species. 


and sub-modules could facilitate quick communication 
by minimizing the local and global energy expenditure in 
communication within the network. The increase in the 
value of fractal dimensions di and fi given by equations 
(1) and (2) (Fig. 3 and Fig. 4) indicates the increase in 
complexity of the network [14, 26]. Since the values of di 


and fi are minimum in monkey brain network as compare 
to the other two species, the modules and sub-modules in 
this network are more ordered and self-organized locally 
as well as globally as compared to the brain networks of 
the other two species. This efficient self-organization in 
the brain network of a certain species might reflect to the 






























FIG. 3: Scaling behavior of modules and sub-modules at vari¬ 
ous levels of organization of the three species, C. elegans, cat, 
and monkey, by calculating network mass (number of nodes) 
as a function of diameter. The right-hand panel is the fractal 
dimension of the three species. 


FIG. 4: Fractal nature of modules and sub-modules at var¬ 
ious levels of organization of the three species, C. elegans, 
cat, and monkey, by calculating the number of intra-edges as 
a function of diameter. The right-hand panel is the fractal 
dimension of the three species. 


fast brain cognition in that species. 

Scaling in centralities and organization 

The betweenness centrality of C. elegans, for the whole 
brain network, modules and sub-modules at different lev¬ 
els, increases as degree of the network increases (Fig. 6 
upper panel) which indicates that hubs in the network 
has significant roles in intra- and inter-modular/sub- 
modular signal processing at different levels. Since high 
value of betweenness centrality of a node of degree k 
reveals that the node could establish quick communi¬ 
cation with other nodes in the network/module/sub¬ 
module through short paths [28-30], hubs in the C. el¬ 
egans brain network may interfere in various network 
regulations and act as a controller of the network. Re¬ 
moving such few hubs emerged in the hierarchical brain 
may cause rewiring of the nodes in the modules and sub- 
modules at various levels that may introduce new hier¬ 
archical topology of modules/sub-modules. The study 
of betweenness centrality of modules and sub-modules 
at different levels Cb as a function of degree k follows 
power-law distribution given by 

CB{e,)-^k^% (4) 

where i = 1,2, 3,4} is the set of power-law exponents 
for different levels indicated by i. The fitted lines on the 


data of modules and sub-modules show nearly parallel 
feature, and the values of ti are in the range [1.23, 2.12]. 
It is also found that Cb of sub-modules increases signif¬ 
icantly with smaller degree k as the level increases (Fig. 
6 uppermost panel (A)). The data of all modules/sub¬ 
modules in a certain level are scaled using one param¬ 
eter scaling method (see Methods section) to a single 
curve and htted with equation (4), and it is found that 
all the four fitted lines on the four different levels are 
approximately parallel (Fig. 6 uppermost panel (B)). 
Similarly, all the data of modules/sub-modules in all the 
levels are again scaled using the same method and fit¬ 
ted with equation (4) (Fig. 6 uppermost panel (C)), and 
the exponent is found to be e = 1.51. This reveals that 
smaller modules have better communication among the 
nodes within each module, and hub/hubs in each module 
has similar roles (as fitted lines with power law on this 
smaller module is nearly parallel with the fitted lines on 
other modules/sub-modules of other levels) with better 
performance (higher value of Cb at smaller value of k). 
However, as the sub-modules reach the motif level (here 
triangular motif), each node in the motif has equal im¬ 
portance (similar hubs due to same degree of each node), 
and therefore Cb of a motif, which is the smallest funda¬ 
mental module, has a single value which is the largest as 
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FIG. 5: Self-similar properties of modules and sub-modules at 
various levels of organization of the three species, C. elegans, 
cat, and monkey, by calculating inter-modular edges of all 
modules and sub-modules as a function of diameter. The 
right-hand panel is the fractal dimension of the three species. 


compared to that of other larger modules/sub-modules 
in the brain network. We did the same process of anal¬ 
ysis to cat and monkey brain network data, and found 
similar behavior in Cb as a function of k given by the 
scaling law of equation (4) (Fig. 6, second and third row 
panels (A), (B), and (C)), and their exponents are found 
to be e = 3.75 and e = 3.02, respectively. The scaling in 
the power-law behavior of Cb indicates fractal behavior 
of the modules/sub-modules at various levels up to the 
motif level. 

Closeness centrality (Cc) is another measure of cen¬ 
trality which describes how quickly an information from 
(by) a node can be propagated (received) to (from) the 
rest of the network, and can be characterized by the in¬ 
verse of average distance between a given node with other 
nodes in the network [31]. The calculated Cc of C. el¬ 
egans as a function of degree k increases as k increases 
(Fig. 6, fourth row panel (A)) which indicates that the 
increase in Cc with k exhibits shorter average path length 
(see equation (9) in Methods) meaning faster information 
processing of the node with the rest of the brain network. 
This means that larger hubs (larger k) are able to com¬ 
municate with the rest of the nodes in the brain network 
of C. elegans faster than the nodes with smaller fc, which 
is true for hubs in modules/sub-modules other than motif 


where every constituting nodes have same k. The data 
of modules and sub-modules at various levels obey the 
following power-law behavior: 

Cc{r]i) ^ ( 5 ) 

where the set {T]i;i = 1 , 2 ,3,4} are closeness central¬ 
ity exponents at various levels. The fitted curves with 
equation (5) on the data of modules and sub-modules at 
various levels of the network are approximately parallel 
(Fig. 6 , fourth row panel (A)). Scaling of modules/sub¬ 
modules at each level are done using one parameter scal¬ 
ing method and then fitted with equation (5) (Fig. 6 , 
fourth row panel (B)), and it is found that the fitted 
lines are approximately parallel. These scaled data of all 
levels are finally scaled with the same scaling method to 
a single curve and the scaled data obey the power law 
given by equation (5) (Fig. 6 , fourth row panel (B)), 
and found ij = 0.16. The same procedure has been used 
to analyze the data of brain network of cat as well as 
monkey. The similar behavior in terms of scaling and 
structural properties are found in both the two species, 
and the two brain data follow the power-law behavior 
given by equation (5) with power-law exponent 77 = 0.28 
and 77 = 0.26, respectively. The results show the fractal 
behavior of closeness centrality which may connect the 
network topology of brain networks to brain functional¬ 
ity. 

Eigenvalue centrality (EC) is in favor of highly corre¬ 
lated nodes (which are usually high degree nodes) with 
rest of nodes in a network, and specific nodes which 
connect central nodes within the network relating to 
global network pattern [32]. EC is characterized by well¬ 
connectedness in a network [33], a smooth enough func¬ 
tion [34] , and is a good measure of spreading (receiving) 
power of information of nodes in (from) the network [35] . 
The calculated EC of the brain network of C. elegans 
(Ce) for the network, modules, and sub-modules at var¬ 
ious levels (see Methods) show increase in its values as 
degree k increases, obeying the following power law, 

CE{S,)--k^\ ( 6 ) 

where {Sf, 7 = 1,2, 3,4} is the set of EC exponents at var¬ 
ious levels. As found in betweenness and closeness cen¬ 
tralities, the fitted lines on the data of C. elegans brain 
network, its modules, and sub-modules are nearly paral¬ 
lel with EC exponents in the range [0.5, l.lj. Similarly, 
it is also found that as one goes towards higher levels, i.e. 
smaller module levels, Cc also increases comparatively. 
We then scaled the data of modules/sub-modules at each 
level, fitted with equation (6), and found that the fitted 
lines on the scaled data are also approximately parallel 
(Eig. 6, seventh row panel (B)) with 5 1.1. We then 

rescaled the data of modules and sub-modules in a single 
one and fitted with equation (6) (Fig. 6, seventh row 
panel (C)) and found 6 0.0.72. 

Similar behavior is found in the brain networks of cat 
(Fig. 6, eighth row panels (A), (B), and (C)) and monkey 
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FIG. 6: Scaling in centrality parameters of brain networks of three species, C. elegans, cat, and monkey: (A) centrality 
measures of all modules and sub-modules at various levels of the brain networks, (B) power-law fits on distribution of the 
centrality measures of each level, and (C) scaled centrality data of all modules and sub-modules into a single curve. The first 
three upper rows of panels are for betweenness centrality, next three rows are for closeness centrality, and last three rows are 
for eigenvalue centrality of the three species. 


(Fig. 6, ninth row panels (A), (B), and (C)) following the 
same scaling power law given by eqnation (6), and found 
the values of 6 to be 0.71 and 1.03, respectively. 

Discussion 

The findings of our study suggest that the fundamental 


working principle of brain (in both lower and higher level 
species) is a system level topological self-organization. 
The fractal nature and scaling properties of these brain 
networks show self-similar organization of various topo¬ 
graphical modules/sub-modules at every levels of consti- 































tution, which may relate to the functional brain organiza¬ 
tion, and energy cost in information transfer within and 
among the levels of organization is minimized. In addi¬ 
tion, the few sparsely distributed hubs are tightly bound 
in their respective module and interfere functionalities 
of their own module, but could not influence rest of the 
modules at various levels in brain networks. In terms of 
inter-modular and intra-modular interaction edges, each 
brain network still show fractal nature which indicates 
systematic self-similar information processing at every 
levels and their interference. The decrease in the values 
of fractal dimension in going from C. elegans to monkey 
(lower to higher species) shows that the organization of 
brain networks (in terms of signal processing, topologi¬ 
cal characteristics, and modular organization) is more or¬ 
dered and self-organized systematically in higher species. 
Such topological properties in brain networks allow effi¬ 
cient information processing, constitution of fractal laws 
in the organization, and controlled behavior of hubs in 
the global network properties. 

The centrality measures (betweenness, closeness, and 
eigenvalue centralities) of brain networks, its modules, 
and sub-modules show increase in their values with de¬ 
gree showing that hubs behave as most influencing nodes 
in the modules/sub-modules they are embedded. These 
hubs act as central in the local module/sub-module, and 
they become local quick information spreader and re¬ 
ceiver in the network. However, removing one hub in 
such situation does not cause the network breakdown 
because of the system level organization of the network 
through modules and smaller modules which are com¬ 
pact with their own fundamental rules. The centrality 
data of the brain network, modules, and sub-modules 
at different levels can be scaled into a single power-law 
behavior showing fractal nature. This exhibited fractal 
nature in the brain network could be the consequence 
of the emergence of a few most influencing hubs in each 
module/sub-module at any level of the network except 
at the level of motif where all the nodes in it have equal 
degree. Therefore, in the brain network, modules, and 
smaller modules, most popular node/nodes always exist 
and they take maximum responsibility in regulating the 
network/module/sub-module. However, these hubs’ in¬ 
terference in the network is controlled (due to limited 
number of links to the modules/sub-modules) in such 
a way that they cannot control the other modules/sub¬ 
modules but can regulate them. 

The scaling properties in brain networks reveal com¬ 
plicated self-organization of the network at various topo¬ 
logical levels, and it could probably explain systematic 
organization of functional modules via weak interaction 
among them. This topographic organization may induce 
the origin of brain functionalities even at the absence of 
few hubs or modules. However, the properties of this 
static network do not fully explain the working princi¬ 
ple of the complicated brain network, its dynamics, and 
functional relationships. The studies on dynamics and 
multi-scaled network approach may highlight further in¬ 


teresting insights on brain organization/reorganization. 

Methods 

Data sources 

In this paper, the connection matrices of (1) C. elegans 
neuronal system, (2) 52 cortical areas in cat species, and 
(3) 71 cortical areas in Macaque monkey species are stud¬ 
ied. 

The C. elegans neuronal connectivity data set is 
adapted from Achacoso & Yamamoto [36], the compila¬ 
tion of which is based on the work of White et al. [37] in 
which the neuronal connection were traced with electron 
microscope reconstructions. Further modifications are 
the removal of 20 neurons in the pharyngeal nervous sys¬ 
tem which have no internal connection information [36] 
and the additional removal of three other neurons (AIBL, 
AIYL, and SMDVL) considering their lack of spatial in¬ 
formation [38, 39]. Finally, 277 neurons sharing 2102 
synaptic connections are considered for further topologi¬ 
cal analysis (data set available at http://www.biological- 
networks.org). 

The cat connection matrix used in this study is de¬ 
rived from the original article by Scannell et al. [40]. In 
their paper, they collected information on the thalamo- 
cortico-cortical connections from many published studies, 
and applied the methods of non-metric multi-dimensional 
scaling, optimal set analysis, and non-parametric cluster 
analysis to derive the connection matrix of the 53 corti¬ 
cal areas and 42 thalamic nuclei. Their connection ma¬ 
trix is relatively weighted (0, 1, 2, and 3) according to 
the connection strength (absent/unreported, weak, in¬ 
termediate, and strong, respectively) between each re¬ 
gion. In this paper, only those connections among the 
52 cortico-cortical areas are studied (after ‘Hipp’ area 
is omitted) [41]. The relative weighting is discarded and 
only the presence or absence of connection is considered 
in the respective adjacency matrix. The resulting final 
matrix has 52 cortical areas and 820 cortico-cortical con¬ 
nections. 

Collecting information from the neuroanatomical stud¬ 
ies, Young [42] applied the method of optimization anal¬ 
ysis to map the cortico-cortical connections between 73 
cortical areas of interest in the entire cerebral cortex of 
Macaque monkey. The connection matrix of Macaque 
monkey used in this paper is also based on the study of 
Young [42], with a modification as mentioned in Sporns 
& Zwi [41] in which two areas of interest (Hipp and 
Amyg) are removed resulting to a total of 71 cortical 
areas with 746 interconnections (data set available at 
https: //sites.google.com/site/bctnet / datasets). 

Graph construction and network parameters 

The connection matrices (adjacency matrices) from the 
above-mentioned data sets are used to generate undi¬ 
rected graphs by using igraph R package [43]. For iden¬ 
tifying communities in these graphs, the leading eigen¬ 
vector spectral graph partitioning method (for which al- 
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gorithm is available in igraph package) is implemented 
[44]. In this method, the modularity term is expressed 
in terms of eigenvalues and eigenvectors of a modularity 
matrix, and the partitioning is done using multiple lead¬ 
ing eigenvectors that optimizes the modularity [45] . The 
communities are then grouped into each topological level. 
For each graph/subgraph we use the Network Analyzer 
[46, 47] and CytoNCA [48] plug-ins in Cytoscape for find¬ 
ing required network parameters such as degree, cluster¬ 
ing coefficient, neighborhood connectivity, betweenness 
centrality, closeness centrality, and eigenvector central¬ 
ity. 

1. Degree distribution 

The degree represents a centrality measure that indicates 
the number of communications a node maintains with 
other nodes in a graph. The degree distribution {P{k)) 
which is the probability that a randomly chosen node 
has a degree k represents an important parameter that 
helps us to identify whether a graph is random, scale free, 
hierarchical, etc. 


m = ( 7 ) 

where Uk represents the number of nodes with degree k 
and N is the total number of nodes in the graph. 


in the network. Thus, for a node i, 


Cc 


n 


( 9 ) 


where d^j represents the geodesic path length from nodes 
i to j, and n is the total number of vertices in the graph 
reachable from node i. 


5. Betweenness centrality 

Betweenness centrality of a node is the measure of the 
extent to which the node has control over the communi¬ 
cation of other nodes. Betweenness centrality (Cb) of a 
node V is computed as follows [50-52]: 

CB(i)= (10) 

where N is the set of nodes, s and t are nodes in the 
graph different from i, is the number of shortest path 
from s to t, and path through i in the case of <Jst(i)- The 
betweenness centrality value is normalized by dividing 
with the number of node pairs (excluding node i). 

6. Eigenvector centrality 

Eigenvector centrality of a node i (vi) in a network is 
proportional to the sum of i’s neighbor centralities [53], 
and it is given by 


£. Neighborhood connectivity 

Neighborhood connectivity of a node i represents the 
average connectivities (average degrees) of the nearest 
neighbors of node i [49]. 

3. Clustering co-efficient 

Clustering co-efficient is a measure of how strongly a 
node’s neighborhoods are interconnected. Graph theo¬ 
retically clustering coefficient is the ratio of the number 
of triangular motifs a node has with its nearest neighbor 
to the maximum possible number of such motifs. For 
an undirected graph, clustering coefficient (Ci) of the ith 
node can mathematicalfy be expressed as 


2ei 

ki{ki - 1 )’ 


( 8 ) 


where Ci is the number of connected pairs of nearest- 
neighbor of the Ah node, and ki is the degree of the ith 
node. 


'^i = \ 

j=nn(i) 

where nn(i) indicates nearest neighbors of node i in the 
network. A is the eigenvafue of the eigenvector Vi given 

by 

Avi = Xvi, ( 12 ) 

where A is the adjacency matrix of the network. The 
principal eigenvector of matrix A, which corresponds 
to maximum eigenvalue Amax, is taken to have positive 
eigenvector centrality scores [35]. 

7. Modularity 

Finally, modularity is the measure of how well a network 
is divided in communities [54]. Modularity {Q) is express 
as follows: 

= (13) 

ij 


Centrality measurement 

In addition, other important centrality measures include 
(1) closeness centrality, (2) betweenness centrality, and 
(3) eigenvector centrality. Centrality measures are help¬ 
ful in identifying influential node(s) in a graph. 

4- Closeness centrality 

Closeness centrality (Cc) of a node is the reciprocal of the 
mean geodesic distance between the node and all other 
nodes reachable from it [31]. Therefore, it represents how 
fast information is spread from the node to other nodes 


where m is the total number of edges in the community, 
Aij is the adjacency matrix of size i x j, k represents 
degrees, and the S function yields I if nodes i and j are 
in the same community. 

Scaling nature of topological parameters 
The data of topological parameters (probability of de¬ 
gree distribution, clustering co-efficient, and neighbor¬ 
hood connectivity) and centrafity parameters (between¬ 
ness, closeness, and eigenvector centralities) of the net¬ 
work, modules, and sub-modules at various levels (Fig. 
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2 and Fig. 6 ) in log-log plot show approximately paral¬ 
lel power-law fit lines. We follow one parameter scaling 
theory [23-25] to scale the data given by 


A 



(14) 


where F’ is a scaling function. For topological pa¬ 
rameters A(fc) = P{k),C{k),Cn{k) and for centralities 
A{k) = C'b(A:), C'c(fc), C'E(fc) with corresponding (p val¬ 
ues after fit. The calculated f after fitting each data 
of network/module/sub-module corresponds to the min¬ 
imum path length of the network/module/sub-module 


approximately. This fitting procedure gives ns F Ki 
constant. Hence, we found the following scaling law: 

A(/c) ss constant x (15) 

where cp = {—a, —/3, — 7 } for C{k), Cn{k), and P{k), re¬ 
spectively, and (p = {e, 7 , <5} for Csik), Cc{k), and C'E(fc), 
respectively. 
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